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Abstract 



We reanalyze the allowed parameters for inhomogeneous big bang nucleosynthesis in light of 
the WMAP constraints on the baryon-to-photon ratio r] and a recent measurement which has set 
the neutron lifetime to be 878.5 it 0.7 it 0.3 seconds. For a set baryon-to-photon ratio r/ the 
new lifetime reduces the mass fraction of ^He by 0.0015 but does not significantly change the 
abundances of other isotopes. This enlarges the region of concordance between ^He and deuterium 
in the parameter space of r] and the IBBN distance scale rj. The ''Li abundance can be brought 
into concordance with observed ^He and deuterium abundances by using depletion factors as high 
as 9.3. The WMAP constraints, however, severely limit the allowed comoving (T = 100 GK) 
inhomogeneity distance scale to rj ~ (1.3 — 2.6) x 10^ cm. 
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I. INTRODUCTION 



Big bang nucleosynthesis (BBN) plays a crucial role in constraining our views on the 
universe. It is essentially the only probe of physics in the early radiation dominated epoch 
during the interval from ~ 1 — 10^ sec. Moreover, big bang cosmology is currently undergoing 
rapid evolution based upon high precision determinations of cosmological parameters and 
improved input physics. Thus, it is important to scrutinize all possible variants of BBN and 
understand the constraints which modern observations place upon their parameters. This 
paper summarizes the current status on one such important variant of the big bang, namely 
one in which the baryons are spatially inhomogeneously distributed during the epoch of 
nucleosynthesis. We consider the constraints from WMAP and the latest observed elemen- 
tal abundances and include recent measurements of thermonuclear reaction rates and the 
neutron lifetime. We show, that although the parameter space for inhomogeneous big bang 
nucleosynthesis (IBBN) is significantly limited, interesting regions remain and this paradigm 
continues to be a viable possibility for the early universe. 



II. BACKGROUND 



At a high temperature in the early universe, e.g. T ~ 100 GK (the corresponding age of 
the universe is ~ 0.01 seconds) baryons mainly exist in the form of free neutrons and protons. 
During this epoch neutrons and protons are rapidly interconverted via the following weak 
reactions m 



n + ^ p + e~ 
n + e"*" ^ p + t'e 

n ^ p + e~ + Ue- 

As long as these reactions are in thermal equilibrium the ratio of neutrons to protons 
is given by (n/p) = exp[—Am/kT] where Am is the mass difference between neutrons and 
protons. When the universe cools to a temperature T = 13 GK these weak reactions fall 
out of equilibrium. The n/p ratio after that time decreases only due to neutron decay. 

When the universe cools to a temperature T «i 0.9 GK the nuclear reaction n + p d 
+ 7 falls out of nuclear statistical equilibrium. Deuterium production becomes significant, 
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leading to the synthesis of increasingly heavier nuclei through a network of nuclear reactions 
and weak decays. When nearly all free neutrons are incorporated into ^He nuclei, ^He 
production virtually stops. The final mass fraction can be approximated as 



Y, ^ 2(n/p)/[(n/p) + 1]. (1) 

Detailed final abundances depend most sensitively on the baryon-to-photon ratio rj. Ob- 
servational measurements of ^He, deuterium and ^Li abundances can be compared with 
abundances calculated from BBN to put constraints on rj. Constraints can also be put 
on 7] from measurements of the Cosmic Microwave Background (CMB). The physics of 
accoustic oscillations in the CMB angular power spectrum has a dependence on various cos- 
mological parameters (See Melchiorri for comprehensive references to methods of CMB 
calculations), including the baryon density factor VL^h?, and therefore 77. Recent CMB mea- 
surements using WMAP Q set fi^/i^ = 0.0224 ± 0.0009 and r/ = 6.13 ± 0.25 x 10"^°. See 
Lara ^] for references to previous attempts to constrain rj using BBN or the CMB. 

The final n/p ratio depends on the neutron lifetime, Xn, which partially determines when 
weak reactions fall out of equilibrium. The most recent lifetime world average |7[ is Tn = 
885.7 ± 0.8 seconds. But a recent measurement by Serebrov et al. ?8| sets = 878.5. ± 
0.7stat ± O.Ssys seconds, which differs from the world average by six standard deviations. 
Serebrov et al. measured the storage loss rate Tg^orage = '^n^ + ''loss °f ultracold neutrons in a 
gravitational trap, and then extrapolated the data down to r^^l^ ^ to get a value for 
alone. Their storage loss rate was as much as a factor of two smaller than measurements 
from previous experiments. So their measurements were much closer to the true neutron 
decay rate, greatly reducing the sytematic error in the measurements. Also, Serebrov et al. 
used the neutron lifetime as a unitarity test of the Cabibbo Kobayashi Maskawa (CKM) 
matrix. The Standard Model of particle physics predicts that the CKM matrix is unitary 
0, Q| . The neutron lifetime by Serebrov et al. was in better agreement with this unitarity 
prediction than the world average lifetime, as shown in Figure 4 of Serebrov et al. Isl. 

n n 

Mathews et al. [HI] applied the measurement Yp = 0.2452±0.0015 from Izotov et al. |12| . 



which is derived from observations of extraga^ 
galaxies, to a Standard BBN (SBBN) code |i 



actic Hll regions in low metallicitv irregular 
13j available for public download '14* . In the 



results using the world average neutron lifetime, the above la range of Yp corresponds to a 



1. Cl 



range r] = 4.8 ± 0.8 x 10 not in agreement with the WMAP results. Mathews et ah 
p.U„ Se.b.ov e. aLMe.i.e v.ue i,.o ..e SBBN code Q^. The new hfetime has the 
effect of lowering Yphj ^ 0.0015 for all values of baryon-to-photon ratio t]. The la range of 



rj corresponding to 1^ = 0.2452 ± 0.0015 1^ shifts to rj = 5.5 ± 0.9 x 10^^°. The new range 
of rj is in agreement with the range rj = 6.13 ± 0.25 x 10^^° from the WMAP measurements. 

This article will discuss how the measured isotope abundances and the WMAP data 
constrain parameters from an IBBN cosmology when using the Serebrov et al. neutron 
lifetime. Various theories of cosmic phase transitions may lead to the formation of bar yon 
inhomqgeneous regions, for example during a first order quark hadron phase transition [if 
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23| or during the electroweak phase transition 



Inhomogeneities can have planar, cylindrical or spherical symmetry 
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barron inhomogeneous regions can be model 
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3 43, Q, 




The inhomogeneous region is broken up into 64 zones, using a stretching 



function by Kurki-Suonio and Matzner j42]. Within each zone s the time evolution of the 



number density n{i, s) of isotope species i obeys the following equation 



dn{i, s) 
dt 



j,k,l 

a Id 

-3-?T,(z, s) H — 

a rP or 



+ NMl ^^^^ 



r^D, 



dn{i, s] 



(2) 



' dr 

The first two terms correspond to nuclear reactions and weak decays that create or destroy 
isotope i within zone s. The weak reactions that interconvert neutrons and protons go 
in these terms for instance. The third term coresponds to the expansion of the universe. 
The final term corresponds to diffusion of isotope i between zones. Neutron diffusion is the 
most prominent diffusion because of the neutrons' lack of charge. The influence of neutron 
diffusion on IBBN results is described in detail by Lara 0]. In both Lara ^ and this article 
the authors use expressions derived by Jedamzik and Rehm 55l 6lL 1621 . |63| for the diffusion 
coefficients. Thermonuclear reaction rates from Angulo et al. j6^J were used in runs for both 
Lara |6] and this article. 
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III. RESULTS 



The following results are for a cylindrically symmetric IBBN model with an initial high 
baryon density in the outermost zones of the model. Parameters used to define the model 
include Rp, the initial ratio of the high to low baryon densities, and /t,, the volume fraction 
of the high density region. Rp is set to 10^, and is defined such that the radius of the low 
density cylindrical core equals 0.925 of the total radius (the high density outer shell then 
has a thickness equal to 0.075). This model was used by Orito et al. |;29j and by Lara 6, 6^ 
to optimally satisfy the light element abundance constraints and also the inhomogeneous 
geometry which is calculated in an effective kinetic model of QCD Final abundance 
results are presented as contour maps which are dependent on rj and the inhomogeneity 
distance scale rj, the size of the outermost zone at the initial temperature of 100 GK. For 
the smallest values of r^, neutron diffusion homogenizes both neutrons and protons early 
enough for the resulting abundances to have the same values as abundances from SBBN 
models. 

Note that the code used by Mathews et al. [l^ uses Bessel function expansions to 
solve for the rates of the neutron proton interconversion reactions 66|. The code used in 
this article uses a numerical integration method to solve for those rates 67[ . The difference 
in the methods of calculating the rates can lead to discrepencies in the value of Yp on the 
order of 0.001. There are also differences in the calculation of the Coulomb correction and 
the zero temperature radiative correction between the codes that affect the results 0, Q ■ 
These differences will be discussed and resolved in a subsequent article. 

Figure n] shows both observational constraints for deuterium in an rj vs. contour map 
and the CMB constraints rj = 6.13 ± 0.25 from the WMAP results P]. The deuterium 
observational constraints D/H = 2.78^Q 3g x 10^^ are the average of measurements of six 
absorption-line systems towards five Quasi Stellar Objects [g^. In the SBBN model the 
deuterium constraints correspond to r/ = (5.6 — 6.7) x 10~^°, encompassing the WMAP 
constraints. As can be seen in Figure Q the results for the smallest values of are equivalent 
to SBBN results. The bends in the deuterium contour lines are explained in detail in Lara . 
In models with from 25000 cm to 10^ cm, neutron diffusion occurs between weak freeze- 
out and nucleosynthesis. Contour lines are shifted to lower r] (lower deuterium production) 
because nucleosynthesis is concentrated in regions with large proton density. For r, from 
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^ 10^ cm to 2.0 X 10^ cm neutron back diffusion does not reach all regions of the model 
during the time of nucleosynthesis, and countour lines shift instead to higher rj (greater 
deuterium production). At ?5i 2.0 x 10^ cm neutron diffusion peaks at the same time 
as nucleosynthesis. The deuterium contour lines shift to lower rj and the IBBN model is 
approximately the average of a high baryon density SBBN model and a low density SBBN 
model. Due to these bends the deuterium constraints are in concordance with the WMAP 
results for three ranges of distance scale: r, < 9000 cm which includes the SBBN case, 
Ti = (1.3 - 2.6) X 10^ cm and = (2.1 - 2.9) x 10^ cm. 

Figure El shows observational constraints for ^He for the cases of both the world average 
neutron lifetime and the neutron lifetime by Serebrov et al., along with the deuterium 
and WMAP constraints. Figure |21 shows Yp = 0.246, the 2a maximum constraint of the 
measurement 0.242 ± 0.002 by Izotov and Thuan l70|l. This measurement is more recent 

n 

than the previously adopted measurement by Izotov et al. [lj| . If the world average neutron 
lifetime is used Yp = 0.246 corresponds to ?7 = 6.1 x 10^^° for the cases of SBBN and IBBN 
with the smallest values of r^. Only a very narrow range 1] = (5.88 — 6.1) x 10^^'^ is then 
permitted for < 2000 cm. The region of concordance for = (1.3 — 2.6) x 10^ cm however 
is contained within this ^He constraint. Note though that the extent of systematic errors in 
^He observations are controversial. Olive and Skillman [^| adopt a conservative approach 
to accessing uncertainties and report a very large range 0.232 <Yp< 0.258 due to correlated 
errors. This range should eventually come down as the number of systems with uncorrelated 
systematic errors increase. 

A shorter neutron lifetime means that more neutrons are converted to protons before 
the onset of nucleosynthesis, and Yp for a set value of t] is smaller. But the lower neutron 
abundance has no perceptible effect (< 1 %) in the final abundances of all other isotopes 11 1. 
Since the neutron lifetime is independent of baryon inhomogeneity, the IBBN run shows the 
same results. For a set r], Yp is reduced by 0.0015 while D/H is unchanged, regardless of 
the value of or r^. The contour line of the ^He maximum constraint 3] shifts to higher rj 
by 1 xlO^^°. Using the world average neutron lifetime, the region of concordance between 
the deuterium and ^He constraints excludes a range (4.5 x 10'^ — 1.3 x 10^) cm of r^. The 
shift in the ^He contour line due to the Serebrov et al. lifetime allows for a concordance 
between deuterium and '^He for all values of < 6.4 x 10^ cm. The WMAP constraints 
are in concordance with the deuterium constraints and ^He constraints with the Serebrov 
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et al. lifetime for ri < 5500 cm. The region of concordance for = (1.3 — 2.6) x 10^ 
cm is still contained within the new ^He constraint. For concordance with the region with 
Tj = (2.1 — 2.9) X 10'' cm a maximum value Yp = 0.259 is needed. This value is close to 
the maximum limit reported by Olive and Skillman, but that maximum limit is expected to 
diminish as the errors improve. 

The primordial abundance of ''Li is in dispute because of disagreement over the calculation 
of the effective temperatures of metal poor halo stars needed to determine the abundances. 
Ryan et al. [t^I and Melendez & Ramirez [3| calculate different values of the effective 
temperature and different ranges of the ^Li abundance. Figure IHl shows the Ryan et al. ^Li 
constraints ^Li/H = 1.23+°:| X 10 along with the deuterium, WMAP and ''He constraints 
when using the Serebrov et al. lifetime. Figure 0] does the same but instead shows the 
Melendez & Ramirez constraints ^Li/H = 2.34!![q gg x 10~^°. Regions of concordance between 
the ^He, deuterium, and WMAP constraints are highlighted in yellow (color online). 

No region of concordance exists between the ^Li constraints of Ryan et al. T^l and the 
constraints from deuterium and ^He for either the SBBN model or the IBBN model. A 
depletion factor due to stellar processes is needed for concordance. A factor of 2.5 would 
bring the ''Li constraints into concordance with the VFM^P constraints r/ = 6.13±0.25xl0^^° 
for the SBBN case and for the IBBN case with distance scale < 2500 cm. But for rj «i 
25000 cm nucleosynthesis in high proton density regions leads to considerably increased 
production of ^Be [6j. The ^Li contour lines shift to lower t] than the contour lines for 
deuterium and ^He. Larger depletion factors can then be permitted in the IBBN model. 
Depletion factors up to 5.9 are permitted using the world average neutron lifetime j^. The 
main effect of using the Serebrov et al. lifetime is to shift the ^He contour lines to higher t]. 
This shift increases the maximum allowed range for the depletion factor, to 9.3 as shown in 
Figure El The constraints by Ryan et al. with this depletion factor will be in concordance 
with the narrow region = (1.3 — 2.6) x 10^ cm permitted by the deuterium, ^He, and 
WMAP constraints. 

Figure m shows a thin region of concordance rj = (5.88 — 6.0) x 10~'° between the WMAP 
constraints and the 2a limit of the ''Li constraints by Melendez & Ramirez [tsI . A depletion 
factor of 1.2 would bring these ^Li constraints into concordance with the whole range rj = 
6.13 ± 0.25 X 10-1° of the WMAP constraints for SBBN and for IBBN with < 3000 
cm. Using the world average neutron lifetime a maximum depletion factor of 2.8 of the 
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Melendez & Ramirez constraints is required by this IBBN model [6|. Using the Serebrov et 
al. neutron hfetime, a maximum depletion factor of 4.47 brings the Melendez & Ramirez 
constraints with the region = (1.3 — 2.6) x 10^ cm permitted by the deuterium, '^He, and 
WMAP constraints. 

The results in this article are for a cylindrically symmetric IBBN model with an initial 
low density core and high density outer region. Changing the symmetry of the model or 
changing the values of parameters Rp and fy have the effect of lessening or amplifying the 
turns in the contour lines, but the characteristics of the contour lines remain the same 
qualitively. Contour lines shift to lower rj and then to higher rj and then back to lower t] as 
the size of the model increases. The particular model of this article is thus representative of 
IBBN models in general. 

We note that new measurements of the ^Li primordial abundance derived from the ratio 
(^Li/^Li) measured in the interstellar medium 3, 13 support a larger range of permit- 
ted depletion factors, consistent with those required here, leaving an enhanced ^Be or ^^B 



abundance as a possible observable signature for the IBBN 



IV. CONCLUSIONS 



We have reanalyzed IBBN in the context of the latest constraints on primordial abun- 
dances, input physics, and cosmological parameters. We have shown that some possible 
IBBN solutions are possible. For example, a cylindrical fluctuation with a comoving radius 
of ~ 10^ cm fixed at a temperature of 100 GK is still consistent with the D/H, Yp and 
WMAP constraints. Depletion factors of up to 9.3 or 4.47, depending on which measure- 
ment is used, can bring the ''Li results into concordance as well. Also any fluctuation up to 
a radius of ~ 10^ cm is viable. Hence, although the standard homogeneous BBN limit is 
still a consistent and the simplest solution, IBBN remains a viable possibility should it ever 
be established that such fluctuations arise from earlier cosmic phase transitions. 
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FIG. 1: Contours consistent with the ±2a measured abundance ratio D/H of deuterium 22| are 







shown in long dashed hues. The CMB constraints from the WMAP measurements laj are also 
shown as solid lines. The WMAP constraint limits the allowed parameter space to three regions. 
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FIG. 2: Constraints from deuterium and CMB observations are shown with the +2a maximum 
constraint of the '^He mass fraction Yp by Izotov and Thuan [,7Q|. The contour for the +2cj maximum 
constraint of Yp is shown by dashed Hne for the world average neutron hfetime [3] and sohd hne 
for the Serebrov et al. hfetime _8]. The region to the right of the maximum constraint is excluded. 
The +2(7 minimum constraint corresponds to values of i] so low it is not necessary to show in this 
figure. The lifetime effect is to shift the Yp contour line to a higher 77 by ~ 1 xlO~^^. 
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:^IG. 3: ( Color Online ) Contours for the ±2(7 abundance constraints on '^Li/H by Ryan et al. 

7^ are shown along with the constraints on deuterium CMB constraint from WMAP j^, and 

yn 
for the case of Tn = 878.5. ± O.Tgtat i O.Sgys ^8]. The region 

of concordance between WMAP, ^He and deuterium is shown in yellow. The contour lines for ^Li 

allow for depletion factors from 2.5 to 9.3 to bring all three observational constraints and CMB 

constraints in concordance with each other. 
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FIG. 4: ( Color Online ) Same as in Figure 01 but with ±20" constraints on ''Li/H by Melendez 
& Ramirez [t^ . A thin region of concordance between the +2a limit of the ^Li constraints and 
the other constraints exists for rj < 1000 cm. A depletion factor of 1.2 improves concordance. 
A depletion factor of 4.47 would bring these ^Li measurements in concordance with the other 
constraints for rj = (1.3 — 2.6) x 10^ cm. 



17 



